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Abstract 

We report progress on a program of gravitational physics experiments using cryo- 
genic torsion pendula undergoing large-amplitude torsion oscillation. This program 
includes tests of the gravitational inverse square law and of the weak equivalence 
principle. Here we describe our ongoing search for inverse-square-law violation at a 
strength down to 10~^ of standard gravity. The low-vibration environment provided 
by the Battelle Gravitation Physics Laboratory (BGPL) is uniquely suited to this 
study. 



1 Introduction 



The torsion pendulum remains the instrument of choice for a variety of labo- 
ratory tests of gravitational phenomena. Such experiments necessarily require 
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the detection of extremely small torques. First employed in the 18*^^ century, 
this instrument was operated in a mode in which the signal torque was detected 
as an angular deflection of a pendulum in the presence of a source mass. This 
"deflection method" is susceptible to certain systematic errors, even a small 
tilt of the apparatus. One way to overcome this drawback is to operate with 
the pendulum executing torsional oscillations and to detect the presence of a 
small torque on the pendulum arising from interaction with the source mass 
through the measured shift in oscillation frequency correlated with changes 
in relative orientation betw een the pendulum equilibr ium position and the 
direction to the source mass! Moore and Bovnton . 1996t ). 



This "frequency method" was commonly used in the 20**^ century for mea- 
surements of the gravitational constant G, but with milliradian oscillation 
amplitude. We operate at a much la rger ampl i tude, near the value that yields 
a maximum in signal-to-noise ratio! Bovnton . 2000). A downside to the fre- 



quency method has been the temperature sensitivity of a torsion fiber's elastic 
modulus. Variation in fiber temperature will produce systematic and random 
errors in the frequency measurement. An effective remedy is to operate the 
pendulum in a cryogenic environment, providing reduced temperature sensi- 
tivity and an opportunity for improved temperature control. 

Details of our cryogenic torsion pendulum system, which uses this frequency 
method, can be found elsew h ere JNewman et al. , 1997 , 19991 : Newman and Bantell . 



T999I : iBantel and Newmanl . l2000a il. The crvogenic environment has a num- 



ber of benefits besides temperature stability, including low thermal noise 
■^^), high fiber Q (> 10^), slower drift of angular equihbrium position, 
easily achieved high vacuum, and effective magnetic shielding with supercon- 
ducting materials. Concerns of bias due to fiber nonlin earitv and anelasticity 



have b een addressed and appear not to be significant (jBantel and Newman 
2000bl) . 



Our apparatus is located at the Battelle Gravitational Physics Laboratory 
(BGPL) on the Department of Energy Hanford site in the arid lands of east- 
ern Washington State. This facility is several kilometers from public access 
and associated anthropogenic noise, and was developed largely because of its 
low ground-motion background (power density less than 10~^^ m^/Hz at 10 
Hz). The photo in Figure 1 shows the entrance to this underground facil- 
ity developed by the UW/UCI collaboration with the support of the Pacific 
Northwest National Laboratory. 

The siting of BGPL in 1994 was determined by a two-year vibration survey of 
a number of candidate locations in Washington State. A former Nike missile 
launch facility on the DOE-Hanford site was chosen to house our laboratory 
not only because it was significantly quieter than all other locations surveyed, 
but also was found to exhibit an exceptionally quiet vibration environment by 
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Fig. 1. Riley Newman (UCI), Roy Gcphart (PNNL) and Paul Boynton (UW) with 
students and staff in backgound at the entrance to BGPL. 



any standard. This attribute, which is essential to the success of our program, 
follows from three circumstances: 1) A low background of natural seismic 
motion because the exceptionally rigid foundation of BGPL is fairly tightly 
coupled to the massive basalt formations that extend east of the Cascade 
Range. 2) The laboratory is situated underground with virtually no structure 
above grade level; consequently, wind buffeting is minimal. The absence of 
trees is a significant advantage for the same reason. 3) There is virtually no 
anthropogenic contribution to ground motion at this location due to highly 
restricted access to the Hanford site. 

This third item is the most important because human activity commonly over- 
whelms the natural seismic background and fixes a substantially higher level 
of ambient ground vibration than would otherwise be present. This low level 
of anthropogenic noise is primarily responsible for the unique suitability of the 
BGPL for carrying out extremely sensitive experiments on the frontier of grav- 
itation physics. Equally important for BGPL is that restricted access within 
the Hanford site perimeter (6 km distant) assures this low background re- 
mains protected from common anthropogenic vibration sources. For this same 
reason, our instrumentation is fully automated as data cannot be taken until 
personnel vacate BGPL. 
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2 Experimental Approach 



Along with tests of the equivalence principle, it is important to verify the weak- 
field limit of general relativity in a manner that is composition independent; 
namely, through tests of the inverse square force law. An inverse square law 
violation (ISLV) could imply the existence of a superposed non-Newtonian 
force that couples approximately to mass. 

Here we describe a null experiment in which a specially configured torsion pen- 
dulum undergoing large-amplitude oscillations in proximity to a source mass 
also of special form, is potentially able to detect inverse square law viola- 
tions approaching 10~^ of standard gravity at a range around ten centimeters 
— nearly a two-order-of -magnitude improveme nt over the current empirical 
limit at that length scale! Moodv and PailJ . 1993[ ). This improvement is accom- 



plished by an experimental configuration that is only second-order sensitive 
to fabrication errors in pendulum and source mass, and not by a substantial 
reduction in dimensional tolerances. 

There are several ways to characterize a violation of the inverse square law. 
The most common is to add a Yukawa term to the gravitational potential, 
yielding an interaction potential 



V=Vg + VY 

= -^(l + ae-/^), (1) 

where G (at r = oo) is Newton's constant, Mthe mass of the object generating 
the potential, a the strength of the ISLV potential relative to gravity, and A 
a characteristic length scale of the interaction. Figure 2 shows the current ex- 
perimental constraint on a as a function of A, and the projected measurement 
capability of our instrumentation. 
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Fig. 2. Composition-independent, experimental, 2(t upper bound on q as a function 
of A. The dashed curve denotes the design goal of this ISLV experiment. 
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While a Newtonian gravitational potential satisfies Laplace's equation in a re- 
gion with no gravitational sources, a non-Newtonian (ISLV) interaction man- 
ifests itself as a nonzero Laplacian of the potential in a source-free region. 
However, because the Laplacian is a scalar and thus spherically symmetric, 
it cannot produce a torque on our pendulum. We therefore designed our ex- 
periment to detect a nonzero horizontal gradient of the Lapla cian, the low- 



est m ultipole-order indication of ISLV in a torsion experiment (|Moore et al 



1994^ . Because we directly detect the horizontal gradient of the Laplacian, our 



experimental results will be independent of the ISLV model. 

The total interaction energy between an object with density profile p(r) and 
the potential energy(l) is given by 

U = I pVd^r. (2) 



A Taylor expansion of the potential energy in cartesian coordinates about the 
origin yields 



U 



/ p 1 + X— + y— + z— + -x^^ + xy^-^ + ...] Vdxdydz (3) 
J \ ox oy oz 2 ox^ oxoy 



dy dz 2 dx"^ dxdy I 



where the derivatives are evaluated at the origin. These terms can be grouped 
into expressions that have the symmetries of the spherical harmonics. This 
gives a concise, orthonormal basis in which to express the interaction energy. 
The general expansion of the interaction energy may then be written as 

U=J2 VnlmMnlm (4) 
nlm 

where 

M^i^ oc J pr^Yrdh (5) 

and the Vnim^ are the appropriate derivatives of the potential evaluated at 
the origin. Grouping the cartesian terms of Equation (3) into the spherical 
harmonics is nontrivial and nor malizat i on of the moments is also challenging. 
This problem was examined by Moord ( 2000[ ). 

The horizontal gradient of the Laplacian corresponds to the multipole term 



311 



^311^311 

d 



oc ^ (VV) j pr^Y^dxdydz , 



(6) 
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We call M311 and the 311 mass and field moments respectively. The 311 
mass moment is manifestly non- Newtonian because the field moment to which 
it couples, V311 , is a derivative of the Laplacian and would be identically zero 
if the interaction were strictly Newtonian. 



There are two primary considerations that drove the pendulum and source 
mass design. First, to take advantage of the multipole expansion, the ratio of 
the size of the pendulum to the distance to the source mass is less than one, 
ensuring the contribution to the interaction energy from each higher order 
is suppressed approximately by that ratio. Second, we purposely eliminated 
the low-order Newtonian mass and field moments of pendulum and source 
mass respectively, up to the order for which the contribution from the nonzero 
moments falls below estimated systematic errors. In particular, since our signal 
has m = 1 symmetry, we null by design all m = 1 Newtonian pendulum mass 
moments from / = 1 to / = 6 and source-mass field moments from / = 1 to / 
= 8. 



This latter constraint was suggested by previous experiments designed to de- 
tect ISLV using torsion pendulum techniques that relied on a pendulum with 
an exaggerated Newtonian moment. Indeed, a high sensitivity to ISLV could 
be achieved using a barbell with large 1—2, m — 2 mass moment. On the other 
hand, this exaggerated Newtonian mass moment would couple to the residual 
I = 2, m = 2 field moments (arising from fabrication errors associated with 
the source mass) to produce a correspondingly large m = 2 systematic effect 
that would mimic a putative ISLV signal in such an experiment. That is, any 
such scheme to increase the ISLV signal moment would necessarily increase 
the limiting systematic effect as well. 



By designing a pendulum that has no low-order Newtonian moments we can 
reduce such systematic effects by two orders of magnitude. Residual Newto- 
nian moments resulting from fabrication errors in the source mass would not 
directly interact with the pendulum since it is specifically designed not to 
detect them. It is therefore only the residual Newtonian mass moments result- 
ing from errors in fabricating the pendulum that can couple to those residual 
Newtonian field moments. As a result, departures from the design configura- 
tion are manifest only in second order in these errors. This design strategy 
dramatically reduces systematic effects arising from standard machining tol- 
erances, metrology, and material density inhomogeneity without resorting to 
heroic fabrication techniques. 
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3 Experimental Design 



The result of pursuing this approach to pendulum design is shown in Figure 
3. This configuration was fabricated at UC Irvine from fused silica parts. It 
has a mass of 240 g and is approximately 8 cm in diameter. The pendulum 
is suspended by a BeCu fiber that, under the weight of the pendulum, is 
stressed to about 60% of its tensile strength. The pendulum is gold coated and 
the fiber electrically grounded to suppress electrostatic torques (important 
because the minimum detectable signal energy for S/N = 1 is only on the 
order of tens of eV). The fabrication tolerance of the pendulum parts and 
their assembly is 6 /im for most dimensions. The parts are bonded together 
using a room-temperature hydroxide catalysis technique that locally dissolves 
the surfaces of the fused silica parts , fusing the m together in a siloxane bond 
and liberating water in the process! Gwol . 1998[ ). These bonds can withstand 
repeated cooling to ~2K when placed in the cryostat. This technique facilitates 
precise component positioning with minimal added mass. 




Fig. 3. Rendition of the ISLV pendulum on left and realization on right (prior to 
coating with gold). 

The ISLV source mass is designed in a fashion similar to the pendulum and was 
assembled at U. Washington. All m = 1 Newtonian derivatives of the potential 
below / = 8 are eliminated. The even I, m = 1 potentials are naturally zero 
because of source-mass inversion symmetry about the horizontal mid-plane. 
Our design therefore has four degrees of freedom used to null the 331, 551, 771, 
and 220 potentials. The 220 potential is eliminated to render the source mass 
insensitive to small tilts about a horizontal axis, which would generate a 221 
moment. The resulting design configuration is shown in Figure 4. The mass is 
assembled by stacking several hundred precision-machined solid stainless steel 
cylinders with a mass totaling about 1500 kg. 

Several trim masses are located at the outer radius of the stack in order to 
remove the residual 221 and 331 field moments resulting from design deviations 
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Fig. 4. Rendition of the ISLV source mass on left and realization on right. 

in the fabrication of the cylinders or their stacking. If necessary, the 441 field 
moment can also be trimmed. The measurement of these residual fields is 
explained below. 

The segmentation of the source mass into these many small cylinders allows 
us to measure and correct for small mass differences resulting from machin- 
ing tolerances and density variations. The ends of the cylinders have shallow 
male/female mating interfaces. As each cylinder must be placed to within 50 
/xm of the design position, a 15 /xm lateral clearance in the interface allows 
for small corrections in the placement of each piece as the stack is assem- 
bled. The precise placement of cylinders is monitored during assembly in real 
time with laser-coordinate-measuring equipment. Hollow stainless steel tubes 
(shown phantom in Figure 4) provide structural spacers and support the cen- 
tral and top portions of the stack. 



4 Experimental Technique 

The torsional motion of the pendulum is tracked using electro-optical tech- 
niques. The source mass, which rests upon an air bearing, is rotated through a 
series of static positions in azimuth. At each position the motion of the pendu- 
lum through several torsional cycles is recorded and fit to a multi-parameter 
model. An azimuthal m = 1 variation in torque on the pendulum would con- 
stitute an ISLV signal and be manifest as a corresponding variation in several 
of these parameters. Commonly measured indicators of a torque would be an 
m — 1 variation in pendulum equilibrium orientation (deflection method) or 
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torsional oscillation frequency (frequency method). Each of these detection 
methods has advantages and drawbacks. 



Prior to the more sensitive cryogenic ISLV search, we are conducting am- 
bient temperature measurements with a second, similar apparatus that re- 
lies on a different scheme to detect ultra-small torques with a torsion pen- 
dulum. We instead measure the m = 1 variation in the amplitude of the 
second ha rmonic of the pendulum motion as our signal ("second-harmonic 
method" ) (|Bovnton . 200Cl( ). This choice is made for two reasons. First, the am- 
plitu de of the second harmonic is insensitive to temperature variations of the 
fiber! Moore and Bovnton . 19961 ). Second, the second harmonic amplitude is 
insensitive as well to tilt of the fiber attachment point. As mentioned ear- 
lier, these effects can be limitations when using the deflection or frequency 
methods. 



The temperature and tilt insensitivities of the second-harmonic method come 
at the possible cost of a longer integration time relative to the frequency 
method if the measurements are dominated by additive white noise (e.g. shot 
noise) rather than thermal noise of the pendulum system. This is partly be- 
cause white noise on the second harmonic amplitude scales inversely as the 
number of periods to the 1/2 and not to the 3/2 power as for the frequency 
method. Additionally, measurement of a signal torque using the second-harmonic 
method is intrinsically nearly an order of magnitude more sensitive to white- 
noise timing error than when the frequency method is used, even for a single 
oscillation period. At room temperature, however, using the frequency method 
would require the flber temperature to be held constant to roughly a fiK to 
achieve the measurement goal of this experiment, a severe challenge. However, 
we plan to repeat the experiment at the Hanford site using the cryogenic ap- 
paratus described earlier. When operating at ~2 K, the reduced temperature 
coefficient of the torsion fiber and improved temperature regulation should 
enable us to reap the strong advantage of shortened integration time. 

The room-temperature pendulum is enclosed in a vacuum chamber evacuated 
to a pressure below 10~^ mbar where the damping of the torsion oscillations 
is dominated by the fiber and not by gas in the chamber. A magnetic damper 
suppresses low frequency (non-torsional) mechanical modes of the pendulum- 
fiber system. A three-layer magnetic shield demonstrably renders magnetic 
effects negligible. 

The residual 221 and 331 mass moments of the pendulum, which result from 
the non-zero dimensional tolerances associated with the fabrication process, 
are suppressed by four movable trim masses attached to the outer face of 
the thicker disk shown in Figure 3. These critical mass moments are first 
measured by stacking the source mass in configurations that exaggerate the 
potentials coupling to these moments. Figures 5(a), 5(b), and 5(c) show the 
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221, 331, and 441 exaggerated configurations respectively. The 441 moment 
of the pendulum is not trimmed, but is measured to verify that it falls below 
the maximum acceptable value. 




Fig. 5. (a) Rendition of the exaggerated 221 source mass, (b) Rendition of the 
exaggerated 331 source mass. The support structure, composed of hollow tubes (in 
phantom) and trays similar to those in the standard ISLV configuration, are shown 
as well as solid source-mass elements, (c) Rendition of the exaggerated 441 source 
mass, (d) Rendition of the exaggerated 221 pendulum. 

The lowest-order field moments of the source mass are also measured and if 
neccessary trimmed to conform to the design goals. The pendulum used for 
this trimming procedure has an exaggerated 221 mass moment and is shown in 
Figure 5(d). This pendulum is used directly to measure and trim the residual 
221 field moment. Moreover, because 

1^31 = ^V22l (7) 

and 

1^441 = ^^221 (8) 

the residual 331 and 441 moments can be measured with the 221 pendulum 
alone by measuring the 221 field at several elevations relative to the mid- 
plane of the source mass. The linear and quadratic variation of the signal with 
elevation then yields the magnitude of the 331 and 441 fields. Once these fields 
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arc measured, they can be trimmed away iteratively using the trim masses 
located at the back of the source mass stack. 



5 The Design Goals 

Uncompensated fabrication errors in the ISLV pendulum and source mass will 
result in spurious apparent values of the measured parameter a in Equation 
1. To assess the resulting error in the measurement of a we first estimate the 
magnitudes of the fabrication errors based on measurements of the completed 
source-mass cylinders, the measured performance of our 3-D laser coordinate 
measuring equipment, and practical experience in manufacturing fused silica 
parts. For the source mass, 50 /im placement errors of the individual cylinders 
dominate their 6 \xm fabrication errors whereas for the pendulum it is the 6 
\xm thickness tolerances of the fused silica components that dominate. 

Departures from designs due to placement and machining errors may be corre- 
lated, so simply combining them in quadrature could underestimate their true 
contribution. Instead, worst-case error estimates for mass and field moments 
were based on these placement and fabrication errors as maximum deviations 
from the three-dimensional design configuration. Once the apparatus is fully 
assembled, and we are able to empirically measure the 221, 331, and 441 mass 
moments and field moments, the measured values should be comfortably less 
than these worst-case estimates. 

To gauge the magnitude of the fundamental constraints posed by systematic 
error from residual Newtonian gravitational coupling between pendulum and 
source mass. Table 1 gives the equivalent contribution to a (at a range A = 12 
cm) for interactions resulting solely from the worst-case fabrication and place- 
ment errors. We see that the hmiting interactions are the 441 and 771 terms 
which yield a contribution to a of a few times 10~^. Aside from these limiting 
systematic errors, there is also the statistical measurement uncertainty. For 
practical integration times at room temperature, we may be able to constrain 
a to less than 10~^. When the experiment is repeated at ~2 K, the thermal 
noise is expected to be significantly reduced, possibly allowing an upper limit 
on a below 10~^. Cryogenic operation may effectively remove noise consid- 
erations and leave the measurement limited by the Newtonian gravitational 
interactions. Such a limit is implied by Table 1, and would yield a constraint 
on a of about 5 x 10~^, hmited by the worst-case 441 interaction. It may be 
possible to trim the 441 moment of the source mass by as much as 90% leaving 
the worst-case 771 interaction as the limiting systematic effect. 

The room-temperature apparatus shown in Figures, 6 and 7 has been assem- 
bled and tested at the University of Washington and is being moved to the 
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Table 1 

Gravitational second harmonic signals arising from the combined fabrication errors 
of the pendulum and source mass. 



Multipole Signal Equivalent a (A = 12 cm) 



221 


8.22 xlO" 


331 


1.32 xlO" 


441 


4.60 xlO- 


551 


3.58 xlO" 


771 


1.31 xlO" 


991 


3.23 xlO" 



BGPL in July 2006 to be reassembled. We expect much improved performance 
in this ultra-quiet site compared to the semi-urban vibration environment of 
our campus laboratory. The Boeing laser metrology equipment will also be 
relocated to enable reassembly of the source mass. This torsion pendulum in- 
strument will be used first to trim and verify the field moments of the ISLV 
source mass by direct gravitational measurements. As described earlier, the 
exaggerated source mass configurations will allow us to trim the pendulum 
mass moments to meet the design requirements also. At this point we will 
attempt an ISLV search with this lower-sensitivity apparatus, followed by in- 
stallation of the ISLV pendulum in the cryostat shown in Figure 8 in order to 
search down to a around 10^^ or below. We also intend t o use the cryog enic 
apparatus for a high-sensitivity equivalence principal test (|Newmanl . 120011 ) . 
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Fig. 6. Graduate student Ricco Bonicalzi adjusting room-temperature torsion pen- 
dulum apparatus in U. Washington lab (magnetic shields absent). 




Fig. 7. Photo showing 221 pendulum (see Figure 5(d)) with vacuum can removed. 
The pendulum is surrounded by mirrors that direct collimatcd laser beams to and 
away from mirror on pendulum. As pendulum rotates, each departing beam sweeps 
across a split photodiode and the time of passage through a fiducial point is mea- 
sured. These "crossing times" are used to determine details of pendulum motion, 
and ultimately the magnitude of any torque acting on the pendulum. 
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